Field-angle-dependent specific heat measurements and gap determination of a heavy 

fermion superconductor URu 2 Si2 
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To identify the superconducting gap structure in URu2Si2 we perform field-angle-dependent spe- 
cific heat measurements for the two principal orientations in addition to field rotations, and theo- 
retical analysis based on microscopic calculations. The Sommerfeld coefficient y(H)'s in the mixed 
state exhibit distinctively different field-dependence. This comes from point nodes and substantial 
Pauli paramagnetic effect of URu2Si2. These two features combined give rise to a consistent picture 
of superconducting properties, including a possible first order transition of H C 2 at low temperatures. 

PACS numbers: 74.70.Tx, 74.25.Bt, 74.25.Op 



It is believed that heavy fermion superconductors 
(SCs), because of their heavy effective mass arising from 
strong electron correlation effects, mostly exhibit uncon- 
ventional pairing states other than a s-wave pairing to 
avoid strong on-site repulsion. Among various known 
materials URu2Si2 is relatively old, discovered in 1985 
and has yet full of mysteries. A phase transition at 
T =17.5 K which was thought to be an antiferromagnetic 
order with a tiny moment [2j , is now under lively debate 
on its origin [3j. Under this so-called "hidden order" 
(HO) the superconducting state appears at T c = 1.3 K. 

The specific heat C{T) cx T 2 at low T HQ and nuclear 
relaxation rate Tf 1 cx T 3 @ suggest a line node gap or 
more accurately that the density of states (DOS) at low 
energy N(E) cx \E\. Up to now, further details of the 
pairing symmetry in URu2Si2 remain unknown: Where 
the line node is if any, or other alternatives. 

Usually by a bulk thermodynamic measurement alone, 
it is impossible to locate the position of nodes on the 
Fermi surface (FS), except for using the field-angle- 
dependent methods 0, & H E, [HI, El. For exam- 
ple, d xy and d x 2_ y 2-w&ve gap structures can be distin- 
guished by rotating an external field H relative to crystal 
axis. The maximum (minimum) of the oscillation ampli- 
tude in the Sommerfeld coefficient j(H) in low-T specific 
heat d,H,[ll| or thermal conductivity k(H)/T [p3 | corre- 
sponds to the anti-nodal (nodal) direction [l3( . Note that 
lim T ^nC(H)/T = j(H) = lTr 2 h 2 N(0, H), where N(Q,H) 
is the zero energy DOS (ZEDOS) at the Fermi level. 

In this Letter, using the field-angle-dependent specific 
heat measurement, we try to clarify the gap structure 
of superconductivity in URu2Si2- For the identification 
of the point nodes located at north and south poles on 
the FS in tetragonal crystal, simple measurement by the 
field rotation is not enough, because the twofold oscil- 



lation pattern due to large tetragonal anisotropy of the 
upper critical field H C 2 hinders the gap structure. There- 
fore, as an extension of the angle-dependent methods, we 
examine the H dependence of j(H) for each field orien- 
tation, as shown in Fig. HI To identify the point nodes, 
we use a specialty of them [8j. As seen in Fig. HJb), j(H) 
behaves differently on H \\ c and H \\ a, because main con- 
tributions to 7(i?) come from the excited quasiparticles 
(QPs) with the Fermi momentum perpendicular to H. 
Namely, J C (H) cx H for H || c where the existing point 
nodes are not sensed by the excited QPs, yielding a slow 
rise in j c (H). This linear lAH) resembles that expected 
for the full gap structure [14|. On the other hand, 'Ya(H) 
for H || a shows a steep increase at lower H, because the 
point nodes are effectively sensed by the excited QPs for 
this H direction. This J a (H) behavior is similar to the 




FIG. 1: ^-dependence of low-T -y(H) for H J c (•) and 
H || a (o), calculated by quasiclassical theory [8J. We plot 
7(ff)/7o vs HjHci for the line node on the equator (a) and 
point nodes on the poles (b), as schematically shown there. 
70 denotes the Sommerfeld coefficient in the normal state. 
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FIG. 2: (Color online) C(H, 9) of URu 2 Si 2 at T = 0.34 K. 
The field-rotational plane is perpendicular to [010] (a) and to 
[001] (b) directions. Inset shows C/T vs. T plot at H=0. 



so-called Volovik y/H dependence [15|. This feature is 
absent in the line node case in Fig. [IJa), because the ex- 
cited QPs feel more or less the nodal structure for both 
H orientations, leading to essentially the same behavior 
la,c(H) oc y/H with slightly different coefficient. 

Through our analysis, it will also become clear that the 
Pauli paramagnetic depairing effect is operative for both 
directions in URU2S12. Previously two experiments point 
out this fact; a strong depression of H c2 [16( and decreas- 
ing tendency of the Maki parameter k 2 upon cooling for 
both directions 



17, 



1 1 81 ] . The following experiment and 
analyses reinforce this, sug gest ing a possible first order 
phase transition at lower T 

Experiment. We have carried out the angle-dependent 
specific heat measurements C(H, 9) on a single crystal 
URU2S12 (100 mg weight), grown by the Czochralski 
method and annealed in high vacuum for 7 days at 950°C. 
The crystal had T c of 1.3 K as shown in the inset of Fig. 
[2j C(H, 9) measurements were carried out by means of 
a semi-adiabatic method in horizontal fields up to 7 T at 
a lower T of 0.34 K 

In Fig. [2] the data C(H, 9) are shown under the field 
rotated in the (010) (Fig. g(a)) and (001) plane (Fig. 
H^b)). The former exhibits a large twofold oscillation, 
which is mainly due to a large anisotropy in H c2 ([IqH^ — 
13 T, /i ff c c 2 ~3 T for T— >0). It should be noticed that 
minima of C(H, 9) in Fig. [2ja) do not correspond to 
nodal directions. Unfortunately, any subtle oscillation 
of C(H, 9) in the ac-plane that could result from nodal 
structures is washed out by the strong anisotropy in H c2 . 
Thus, we have to examine the //-dependence of C(H,9) 
at each field orientation. 

We have also searched for a possible line nodal struc- 
ture running parallel to the c axis, which can be detected 
by rotating H within the (001) plane. As shown in Fig. 
[5Jb) there is no apparent oscillation within the experi- 
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FIG. 3: (Color online) (a) Field variation of C(H) of URu 2 Si 2 
at T = 0.34 K for H\\a (open circles) and H\\c (solid circles). 
Solid lines are guides to the eyes, (b) //-dependent part of C 
as a function of H / H C 2- The dotted line is a linear extrapola- 
tion of the low- field part for H \\ c. (c) H- variation of C in the 
normal state at 1.52 K for H || a and H || c. (d) T-variation of 
M/H in the normal state below 25 K for H\\a and H \\ c. 



mental accuracy. We could conclude that there is no line 
node parallel to the c axis because the expected oscilla- 
tion amplitude (3 ~ 4%) [H, Q is beyond the experimen- 
tal resolution (~ 0.3%). However this cannot exclude the 
possibility for the line node along the equator on the FS 
and weaker gap anisotropy parallel to the c axis. 

In order to investigate the possible nodal structure, we 
examined the field variation of the specific heat in detail. 
Figure [3]Ja) displays the results at T=0.34 K for both 
H || a and H \\ c. C C (H) for H \\ c peaks at B c c2 (-2.5 T) 
and reaches the normal-state value at higher H. For 
H || a, in contrast, C a {H) continues to increase at 7 T be- 
cause H^ 2 (~ 12.5 T) is much higher in this direction. Ac- 
cordingly, a significant anisotropy in C a . c {H) develops at 
high H. Surprisingly, this anisotropy in C a . c (H) rapidly 
diminishes with decreasing H and both curves almost 
fall on top of each other below 0.5 T. This unexpected 
behavior can be seen more clearly in Fig. [3]Jb), where 
the field dependent part of C is plotted as a function of 
H/H c2 for the two directions. It is remarkable to find the 
distinctively different field dependence of C ayC (H / H c2 ) at 
low H. The rapid rise in C a (H j H c2 ) is reminiscent of a 
nodal structure. C C (H/H c2 ) shows, in contrast, a much 
weaker linear rise, suggesting a full gap feature. This 
distinctive C{H) behavior for two directions never oc- 
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curs when the line nodes exist on the FS, as discussed 
in Fig. QJa). Thus, the experimental characteristics of 
URu2Si2 in Fig. [3{b) are qualitatively consistent with 
the gap structure with the point nodes locating at the 
north and south poles on the FS, as shown in Fig. [ljb). 

Since superconductivity in URu2Si2 coexists with HO 
that sets in at 17.5 K, one might suspect that its 
field-direction dependent excitations could result in the 
anisotropic behavior of C^ C {H / H C 2) at low H. In or- 
der to rule out this possibility, we carried out the C(H) 
measurements at 1.52 K just outside the superconduct- 
ing state but still well inside the HO phase (Fig. EDJc)). 
For both field directions, C aiC (H)/T is only weakly in- 
dependent; a slight linear decrease by 2~3% has been 
observed at 2 T. Since the C^ C (H)/T variation in the su- 
perconducting state (Fig. [U[a)) is much larger (40^100% 
at 2 T), this change in the normal-state background is 
negligible. In order to back up this argument, we have 
also measured the magnetization M(T) at /xo-ff=0.5 T 
(Fig. EJd)). It is well known that M a (T) for H \\ a is 
nearly T-independent whereas M C (T) for H || c exhibits a 
kink upon HO at 17.5 K and rapidly decreases at lower 
T0. For both directions, M/H below 10 K, where HO 
is well developed, is virtually T-independent. The slight 
upturn in M ajC (T) seen below 5 K in Fig. OHd) is pre- 
sumably due to impurity effect. Then making use of a 
Maxwell relation d 2 M/dT 2 = T^dC/dH, field varia- 
tion of C (H) is expected to be small, consistent with the 
results in Fig. [3fc) . Thus we conclude that the contribu- 
tion of HO to C(H) is small and negligible; the distinc- 
tively different field dependence of C ajC (ii / H C 2) in Fig. 
EJb) arises from superconductivity in the system. 

Compared with theoretical curves in Fig. [TJb), ^(H) 
of the experimental data in Fig. [3jb) is suppressed at 
the middle-ii region for both directions. Especially for 
H || c, starting from a gradual linear increase (dotted 
line) at low iJ, 'y(H) shows concave curvature at high 
H, which is in contrast to the theoretical curve in Fig. 
[]Jb) showing convex curvature. We may attribute the 
anomalous behavior to the Pauli paramagnetic depairing 
that suppresses superconductivity eminently at high H. 
Then H c2 is reduced, and j(H) rapidly rises at high H. 

According to the M(H) measurements [13, EH, the 
Maki parameter k 2 oc 1/ \j (dM/dH) H=H ^ is an increas- 
ing function of T for both orientations, which is opposite 
to an ordinary SC [20| . At H = H C 2 they show large para- 
magnetic moments build in the mixed state. These facts 
clearly indicate that the paramagnetic effect is an impor- 
tant ingredient in fully understanding this material. 

Theoretical calculation. To support above suggestions, 
we calculate the ZEDOS iV(0, H) by the microscopic 
quasi-classical theory in the clean limit including the 
paramagnetic effect. As is done in Ref. [2l|, we selfcon- 
sistently determine the spatial structure of the pair po- 
tential A(r) and the vector potential A(r), to appropri- 
ately evaluate the vortex core contribution, through the 
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FIG. 4: (Color online) if-dependence ofy(H) for H \\ a in the 
linear point node case 0(k) = sin (a) and the quadratic point 
node case (j>(k) = sin 2 6 (b). Solid circles are the experimental 
data in Fig. [^b). /i =0 (A), 0.85 (B), 1.7 (C) in (a), fi =0 
(A), 0.6 (B), 1.0 (C) in (b) 
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FIG. 5: (Color online) The same as in Fig. |4] but for H \\ c. 
fi =0 (A), 1.6 (B), 2.5 (C) in (a), fi =0 (A), 1.3 (B) in (b). 

quasi-classical Green's functions g(ui n , k, r), /(w ra ,k, r), 
and /'(w„, k, r), which are calculated by the Eilenberger 
equation 

{uj n + UiB + k • (V + iA)} / = A0(k) 5 , 

{cj n + i/iS - k • (V - iA)} /t = A*0* (k) ff , (1) 

with g = (1 — fp) 1 / 2 , Keg > 0, Matsubara frequency 
u> n = (2n + l)irT, and effective Zeeman energy [iB. 
/i characterizes the strength of the paramagnetic effect. 
The self-consistent calculation is performed at T = 0.1T C 
in the triangular vortex lattice for k on an isotropic 
Fermi sphere. The local DOS at an energy E is given by 
N(E,r) = (Reg(iuj n -> E + i0+,k,r)) k with the Fermi 
sphere average (■ • By the spatial average of N(E = 
0, r), we obtain the ZEDOS N(0, H) = (N(E = 0, r)) r . 

In Figs. 4 and 5 we examine two cases of point nodes, 
(i) quadratic point nodes by the pairing function </>(k) = 
sin 2 #, and (ii) linear point nodes by </>(k) = sin 9, where 
9 is the polar angle from the c axis. The former (i) is 
consistent with the experiments C(T) cx T 2 0, 0| and 
T^ 1 cx T 3 [6]. However, as shown in Fig. [5fb), the re- 
sulting ^ C {H) exhibits upward curvature at low H, which 
fails to explain the linear ii-dependence of our experi- 
mental data. The latter (ii), used in the calculation for 
Fig. [U can reproduce ^/H-Mke behavior for H \\ a (Fig. 
2{a)) and ii-linear behavior for H || c (Fig. [Ha)). For 
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better fitting to the experimental data, we have to include 
the paramagnetic effect by fi. The best fittings at low H 
are attained by [M a — 0.85 (p, c — 1.60) for H\\a (H || c). 
The discrepancy at higher H is due to thermal effect be- 
cause the experiment was not done at sufficiently low T. 
Near T C (H), C{H)/T is larger than j(H). According 
to the M(H) measurements 17, [l8j . the paramagnetic 
moments at H C 2 are comparable for two directions, i.e. 
HaHc2 — Mc^c2- This indicates fi c / fi a — 2.1, which is 
consistent to our choice fi c /^a = 1.60/0.85 ~ 1.9. 

Our calculation can show that above a critical value 
/i* ^0.4 the first order transition occurs at H C 2 in lower 
T and higher H for a Fermi sphere and s wave pairing. 
Thus, our assigned fi values are within the first order 
region. Since the C(H) data were taken at T = 0.34 K, 
there is no indication for it. However, recent thermal 
conductivity measurements (l9| clearly shows a jump of 
k(H) at H C 2 for both directions around T ~ 0.1 K. This 
is in accord with our assignment for the fj, values. 

Among the possible pairing function in the group the- 



oretical classifications for a tetragonal crystal [22j, [23( 

— \~ 17~y k y 



spin-triplet symmetries such as r x k x + T y k y , 
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i~ z (k x + iky), and (r x + ir y )(k x + ik y ) have point nodes 
fa = iaiUj), but they are unlikely to be realized in 
URU2S12 because the paramagnetic effect is present in 
both directions. In spin-singlet symmetries, k x + k y (Ai s ) 
has two quadratic point nodes, and k z (k x + ik y ) (E s ) con- 
sists of two linear point nodes and a line node. Our anal- 
ysis supports linear point nodes, and excludes line nodes. 
Thus, if we choose k z (k x + ik y ) as a plausible gap func- 
tion with linear point nodes, the line node contribution 
may be smeared out due to the actual FS topology [24| . 
While we have to consider realistic FS topology for con- 
clusive evaluation of the gap structure, the linear point 
nodes in the polar direction are plausible, because low-iJ 
behavior of "f(H) is governed by the low-E QPs excited 
around nodes [la ]. 

CeCoIns is also considered to be a Pauli limited SC 
with a tetragonal structure. The phase diagrams in H vs 
T for both field directions are similar and exhibit Pauli- 
limited behavior. It is understandable because M values 
at H C 2 are almost same for both H \\ c and H \\ a [25j |. 
which matters the paramagnetic effect. This situation 
is the same as in URri2Si2 as mentioned before. These 
two compounds belong to a clean limit SC free from dirt 
effects, such as spin-orbit scattering that tends to mask 
a first order transition due to the paramagnetic effect. 

In summary, we performed field-angle-dependent spe- 
cific heat measurements in URu2Si2- The Sommerfeld 
coefficient 7 (if) exhibits distinctively different behavior 



for H || c and H \\ a. These reveal, supported by micro- 
scopic calculation, that the gap structure posses point 
nodes at north and south poles on the Fermi surface and 
moreover the Pauli paramagnetic effect is important in 
this system. This also suggests a first order transition 
at H C 2 in lower temperatures, which has been confirmed 
recently (l9j |. 

We thank Y. Matsuda for informative discussions. 
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